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We analyze the nonlinear voltage dependence of 
elelctrochemical capacitance for nano-scale conduc- 
tors. This voltage dependence is due to finite density 
of states of the conductors. We derive an exact ex- 
pression for the electrochemical capacitance-voltage 
curve for a parallel plate system. The result suggests 
a quantum scanning capacitance microscopy at the nano- 
scale: by inverting the capacitance- voltage expression 
one is able to deduce the local spectral function of the 
nano-scale conductor. 
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It has been well known that density of states affects 
the capacitancaj of a system. In this work we investigate 
a nonlinear bias voltage dependence of the electrochem- 
ical capacitance of a model capacitor at the nano-scale 
for which the density of states (DOS) also plays the cru- 
cial role. The physical origin of this bias dependence 
which we have examined, is not because of depletion of 
charges on a capacitor plate which has been understood 
in semiconductor research, but because of the finite DOS 
of the plates. We show that the accumulated charge 
on a conductor has a nonlinear bias dependence due to 
DOS effects, it is this nonlinear charge which leads to 
the nonlincaj- capacitance. For conductors at nano-scale, 
it is knowrn that the very small DOS plays an impor- 
tant role for the behavior of capacitance. However the 
DOS induced nonlinear capacitance of nano-conductors 
has never been studiecu before and our investigation sug- 
gests a very interesting application which is the quantum 
scanning capacitance microscopy. 

Central to the problem is the determination of charge 
Pa accumulated on a conductor labeled by a. In general 
Pa is a nonlinear function of bias voltages {Vg}, 

Pa = 2_^ CapVfj + - 2_^ Cap-^VjiV^ + ... 
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where Cq^^ = dvf,Pa is the _nsual electrochemical 
capacitanceB'Q, Cap-y = dv^dv^Pa^ etc are the nonlinear 
electrochemical capacitance coefficients, and CapUVa}) 
is the general voltage dependent nonlinear capaci- 
tance. So fcuj|-investigations on quantum correction to 
capacitanceEl'crQ have only considered the linear term, but 
in this work we focus on a general nonlinear expression 
for CaMV^}). 



To be specific we consider a model parallel plate capac- 
itor connected to electron reservoirs Ly perfect leadscl. 
We adopt the dynamic point of vieTO to calculate the 
electrochemical capacitance, hence a finite bias applied 
at the reservoirs injects a charge density into the capac- 
itor plates which, through interaction, induces a local 
response. The total net charge (the sum of injected and 
induced charge) at a plate a, pa, is thus established at 
equilibrium. We calculate pa by extending the standard 
nonequilibrium Green's function (NEGF) technique. To 
save space we refer interested readers to literature for 
standard technical details, here we will only discuss the 
most important extension to the standard approach. 

The quantum scattering (injection of electrons) is de- 
termined by the retarded and advanced Green's func- 
tion G^'"'{E,U): note that we have explicitly included 
the electro-static potential build-up inside our capacitor, 
U = [/(r), into, these Green's functions. In the Hartree 
approximatioro. 



G''^''{E,U) = 
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where H is the Hamiltonian for our nano-scale cppduc- 
tors written in the familiar second quantized formcl, S'''" 
is the self-energy describing the coupling between the 
conducjtare to the leads which is calculated in standard 
fashioroEj. At Hartree level we determine U{r) by the 
self-consistent Poisson equation 



V^C/ = Aniq / idE/2Tr)G<{E, U) 



(3) 



where the right hand side is just the total net charge 
distribution in our conductors. Within Hartree approxi- 
mation, 

G<{E,U)^G-Y.'Tp{E~qVp)f{E-qVp)G'' , (4) 



where T p is the voltage dependent coupling jjiarameter 
between probe (3 and the scattering regiorclllO. The 
self-consistent equations (pi WM completely determine 
the nonlinear physics at the Hartree level. We empha- 
size again that the important depaiture of our theory 
from the familiar NEGF analysisEl M is that we explic- 
itly include the internal potential landscape U{y) into the 
Green's functions self- consistently. 

The net charge pile-up on a capacitor plate, measured 
from the equilibrium background, is derived from the 
right hand side of Eq. (P). In the wideband limitEJ 
(Fq = constant) 
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where G^'" are the equihbrium Green's functions. Next, 
we formally expand G*" and G° in a power series of the 
internal potential U , and expand the Fermi function in 
series of the bias voltage Va- Collecting terms accord- 
ing to the powers of U and V^, Eq. d|) reduces to the 
following infinite series which can be exactly summed, 
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where the quantity dua/dE is defined as 
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The physical significance of the quantity dua/dE can be 
identified as the„linear spatial dependent local partial 



density of state 
obtained beforeE3t 



(LPDOS). Expression (0) has been 
For a conductor which is weakly 
coupled to external leads, LPDOS gives the local DOS 
of this conductor. Hence the spectral function (Ja{E) 
characterizes the local electronic structure of a nano-scale 
conductor. 

The nonlinear charge distribution gives the general 
electrochemical capacitance versus voltage curve. 



G - [cj^{E + Vi- Ui) - cTi{E)]/{Vi ~ V2) 



(8) 



where we have set electron charge q to be unity. To 
determine C we must obtain internal potentials Ui and 
U2 at the two plates. Far this purpose we introduce the 
geometrical capacitanceQ Go 



Go = [cTi{E + Vi- Ui) - 'Ji{E)]/{Ui - U2) 



(9) 



and 



Go = -[<J2{E + V2- U2) - a2{E)]/{Ui ~ U2) (10) 

For a parallel plate capacitor, Go = A/(47ra) where A is 
the area of the plates and a is their separation. In gen- 
eral Go can be calculated numerically. The two equations 
(p|,p^ determine the internal potentials Ui and U2 when 
the scattering spectral function <Ta and Go are known. 
Eqs. (pi H, M) gives, for the first time, the general elec- 
trochemical C-V curve for quantum capacitors. It is uni- 
versal in the sense that system specific parameters only 
appear in the scattering spectral functions of the conduc- 
tors. From these general results, two useful applications 
follow: 

C-V curve . Our general results allow us to predict 
capacitance-voltage curves. For this application we note 
that there are well established methods for calculating 



the scattering LPDOsllio: either by Eq. (@) after eval- 
uating the Gceen's functions, or using the scattering 
wavefunctionsEJ -0 doa/dE — \ijj\'^ /hv, where v ^ \fE 
is the velocity and h the Planck constant. Hence by solv- 
ing a quantum scattering problem one obtains a a- Let's 
consider a case where f/) is not very sensitive to £', thus 
daa{E)/dE « ba/{2^/E) or (7a{E) = baVs where ha is 
a constant. For this LPDOS, solving Eqs. (||, §, |l^) we 
obtain 



CiVi - V2) 
where 



l3-^l3^-4{bi-b-i)b'ib-i{V,-V2) 
2{bl~bl){Vi-V2) 



(3 = 47ra6f5^ + 26162(^1 + 62)VS 
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The inset of Fig. (1) shows this electrochemical C-V 
curve as a function of V = (V2 — Vi) for two sets of pa- 
rameters ba . The physical reason that G changes with V 
is because the plates have finite DOS. Indeed, by making 
DOS very large (61,62 — + 00) the voltage dependence of 
( |ll| ) disappears and G becomes purely geometrical. Fur- 
thermore ,_it can be confirmed that forropla (11) recovers 
the lineara and second order nonlineard capacitance co- 
efficients when we take the (V2 — Vi) -^0 limit. 
Quantum scanning capacitance microscopy. Our 
general results suggest a quantum capacitance mi- 
croscopy (QSCM). This idea naturally follows from the 
results presented above: since the electrochemical capac- 
itance varies with bias due to a finite DOS of the conduc- 
tors involved, we should be able to find the DOS by mea- 
suring G. Essentially we wish to obtain spectral function 
(J2{E) or local density of states da2/dE as a function of 
energy for an unknown conductor, from a knaatn ai{E) 
of our QSCM "tip" which has been calibratedEj. As the 
QSCM tip is scanned along the surface of a nano-scale 
conductor, or along the surface of a planar dielectric layer 
with nano-conductors buried underneath, experimentally 
one can measure the G( V) curves at each spatial position. 
From Eq. (pi) , we obtain Ui as a function of potential 
difference V using the known cti and the measured G(V), 
by solving the equation cti(1 — Ui) — cri(l) = VC{V), 
where we have set Vi = 0, V2 = V, and E ^ 1 \s set at 
the Fermi energy of the QSCM "tip" (G(V) is measured 
at Fermi energy of the "tip" E\,, e.g. Eq. (O)). Next, 
From Eq.(g) we obtain U2{V). With Ui and U2 we finally 
find (72 (-E) from Eq. dlQ). In particular we can solve the 
spectral function CT2 by representing it into a polynomial: 
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where the coefficients y„i = ((i™cr2/dV™)/(m!) are just 
the linear (m = 1) LPDOS and nonlinear (m > 1) 
LPDOSO. They are obtained by solving the following set 
of linear algebraic equations which come from Eq. (110) , 



(14) 



m=0 



where j = 1, 2, 



Hence by making experimental 



measurements at n different voltages Vj, we obtain the 
functional form of (T2{E) from dlq). Fig. (1) demon- 
strates the principle of QSCM. We use Eq. ( |ll|) as the 
experimentally measured C-V curve (the inset) to simu- 
late a measurement. Then using the QSCM "measured" 
(72 (£^) from Eq. (13), we plot the local density of states 
da2/dE versus energy E. The solid line is the exact 
da2/dE = &2/(2\/^) and the dots are the QSCM re- 
sult. We used 10 voltages in solving Eq. (|lj) and the 
outcome is quite good, while using 3 voltages it already 
represents a rough trend. 

In summary, we have developed a general nonlinear 
DC theory which is applied to investigate, for the first 
time, the full nonlinear charge distribution in nano-scale 
conductors. We have derived an exact expression of the 
electrochemical capacitance versus external bias voltage 
curve for quantum capacitors. This result is generic in 
the sense that all system specific information are included 
in the scattering local density of states. Hence the C-V 
formula has a wide range of applicability. By invert- 
ing this formula, we propose a novel QSCM. The QSCM 
extends thje-, ability of the usual scanning capacitance 
microscopyll3: QSCM includes the quantum corrections 
to capacitance in mapping out the the spatial charge dis- 
tribution; and it gives the local density of states as a 
function of electron energy. We believe such an idea can 
be readily implemented in a scanning apparatus using 
tiny tips as the calibrated conductor, thus allowing mea- 
surements of electronic properties for other conductors 
at the nano-scale. 
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FIGURE CAPTION 

Fig.(l) Operation of QSCM: comparison of "measured" 
LPDOS to the exact one. Solid lines: exact; solid 
circles: fitted LPDOS using 10 voltages; solid tri- 
angles: fitted LPDOS using 3 voltages. Upper set 
of curves are for 62 = 0.002, lower set 62 = 0.001. 
The QSCM tip has been fixed with 61 = 0.003. 
The energy unit is the Fermi energy E^ of the tip. 
Inset: the electrochemical capacitance versus volt- 
age curve for the two sets of 62: upper curve is for 
62 = 0.002. The unit of V is E\,/e. 
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